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ABSTRACT 

Traditionally runaway stars are O and B type stars with large peculiar velocities. We want 
to extend this definition to young stars (up to « 50 Myr) of any spectral type and identify those 
present in the Hipparcos catalogue applying different selection criteria such as peculiar space 
velocities or peculiar one-dimensional velocities. 

Runaway stars are important to study the evolution of multiple star systems or star clusters as 
well as to identify origins of neutron stars. 

We compile distances, proper motions, spectral types, luminosity classes, V magnitudes and 
B — V colours and utilise evolutionary models from different authors to obtain star ages and 
study a sample of 7663 young Hipparcos stars within 3 kpc from the Sun. Radial velocities 
are obtained from the literature. 

We investigate the distributions of the peculiar spatial velocity, the peculiar radial velocity 
as well as the peculiar tangential velocity and its one-dimensional components and obtain 
runaway star probabilities for each star in the sample. In addition, we look for stars that are 
situated outside any OB association or OB cluster and the Galactic plane as well as stars of 
which the velocity vector points away from the median velocity vector of neighbouring stars 
or the surrounding local OB association/ cluster although the absolute velocity might be small. 
We find a total of 2547 runaway star candidates (with a contamination of normal Population 
I stars of 20 per cent at most). Thus, after subtraction of those 20 per cent, the runaway 
frequency among young stars is about 27 per cent. 
We compile a catalogue of runaway stars which will be available via VizieR. 

Key words: catalogues - stars: early-type - stars: kinematics. 



1 INTRODUCTION 

Almost fifty years ago, Blaauw (1961) found that many O and B 
type stars show large peculiar space velocities (> 40km/s). For 
that reason they were assigned the term "runaway" stars. Many 
studies concerning O and B type runaway stars have been published 
since then covering different selection methods (e.g. Vitrichenko 
et al. 1965; Cruz-Gonzalez et al. 1974; Gies & Bolton 1986; Stone 
1991; Moffat etal. 1998). 
Two theories on the formation of runaway stars are accepted; 

(i) The binary supernova (SN) scenario (BSS) (Blaauw 1961) 
is related to the formation of the high velocity neutron stars; The 
runaway and the neutron star are the products of a SN within a 
binary system. The velocity of the former secondary (the runaway 
star) is comparable to its original orbital velocity. 
Runaway stars formed in the BSS share typical properties such as 
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enhanced helium abundance and a high rotational velocity due to 
mass and angular momentum transfer during binary evolution. The 
kinematic age of the runaway star is smaller than the age of its 
parent association. 

(ii) In the dynamical ejection scenario (DES) (Poveda et al. 
1967) stars are ejected from young, dense clusters via gravitational 
interactions between the stars. 

The (kinematic) age of a DES runaway star should be comparable 
to the age of its parent association since gravitational interactions 
occur most efficiently soon after formation. 

Which scenario is dominating is still under debate; however, both 
are certainly taking place (one example for each identified by 
Hoogerwerf et al. 2001; BSS - PSR B 1929-1-10/ ( Oph; DES - AE 
Aur/ ^ Col/ 1 On). 

The selection criteria for runaway stars of previous studies 
were either based on spatial velocities (e.g. Blaauw 1961), tan- 
gential velocities (e.g. Moffat et al. 1998) or radial velocities (e.g. 
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Cruz-Gonzalez et al. 1974) alone. According to Stone (1979) the 
velocity distribution of early type stars can be explained with the 
existence of two different velocity groups of stars: a low veloc- 
ity group containing normal Population I stars and a high velocity 
group containing the runaway stars (see, e.g.. Fig. 1). Since both 
groups obey a Maxwellian velocity distribution also runaway stars 
with relatively low velocities exist. We want to combine previous 
methods in order not to miss an important star because the radial ve- 
locity may be unknown (hence no spatial velocity) or its tangential 
or radial velocity component may be significantly larger than the 
other component (hence one would miss it in one direction). More- 
over, we also want to identify the lower velocity runaway stars by 
searching for stars that were ejected slowly from their parent clus- 
ter. Furthermore, we want to use the term "runaway" star not only 
for O and B type runaway stars but all young (up to ~ 50 Myr) 
runaway stars to account for the possibility of less massive com- 
panions of massive stars (which soon after formation explode in a 
SN) and low-mass stars in young dense clusters which also may be 
ejected due to gravitational interactions. 

In section 2 we describe the selection procedure of our sample of 
young stars. In section 3 we apply different identification methods 
for runaway stars and assign runaway probabilities to each star. A 
catalogue containing our results will be available via VizieR soon 
after publication. 
We give a summary and draw our conclusions in section 4. 



2 THE SAMPLE OF YOUNG STARS 

We start with all stars from the Hipparcos catalogue (Ferryman 
et al. 1997), 118218 in total, and collect spectral types as well as 
V magnitudes and B ~ V colours from the catalogue. According 
to the errata file provided with the Hipparcos catalogue, we correct 
erroneous spectral types. From the new Hipparcos reduction 
(van Leeuwen 2007) we obtain parallaxes (tt) and proper motions 
(Ha ~ l^ct cos 5, fis). We restrict our star sample to lie within 3 kpc 
from the Sun (tt — cr^ ^ 1/3 mas with a^ being the la error on 
-k)} Furthermore, we remove stars in the regions of the Large and 
Small Magellanic Clouds having accidentally tt — cr^ ^ 1/3 mas 
(cf. Hohle et al. 2010). This gives us an initial set of 103217 stars. 
In the cases where the Hipparcos catalogue does not provide 
sufficient spectral types, we take the spectral type from either the 
Simbad or VizieR databases (Skiff 2009; Kharchenko & Roeser 
2009; Abrahamyan 2007; Kharchenko et al. 2007; Bobylev et al. 
2006; Wright et al. 2003; Buscombe & Foster 1999; Grenier 
et al. 1999; Garmany et al. 1982). Missing B ~ V colours were 
amended from different sources (Simbad; available via VizieR: 
Kharchenko et al. 2007, 2004; Bobylev et al. 2006; Zacharias et al. 
2005; Wright et al. 2003; Egret et al. 1994). 

Having now collected parallaxes, spectral types as well as V 
magnitudes and B — V colours we calculate luminosities (L) and 
obtain effective temperatures (T^ff) from spectral types according 
to Schmidt-Kaler (1982) and Kenyon & Hartmann (1995). The 
extinction Ay is determined from the apparent B — V colour 
and the spectral type.^ The initial sample contains 1721 stars also 



included in the list of Hohle et al. (2010) who used additional 
colours to determine L. Thus, for those stars, we adopt their 
luminosities (without Smith-Eichhorn correction). 
For having full kinematics, we obtain radial velocities from 
Simbad or the following VizieR catalogues: Zwitter et al. (2008), 
Kharchenko et al. (2007, 2004), Bobylev et al. (2006), Gontcharov 
(2006), Malaroda et al. (2006), Karatas et al. (2004), Barbier- 
Brossat & Figon (2000), Grenier et al. (1999), Reid et al. (1997) 
and Evans (1967). 



2.1 Age and mass determination 

Since we look for young stars, our first goal is to obtain star ages. 
For that reason we check our star sample for pre-main sequence 
stars and find 236 of them either in catalogues (Chauvin et al. 
2010; Ducourant et al. 2005; Hernandez et al. 2005; Wichmann 
et al. 2005; Valenti et al. 2003; Kohler et al. 2000; Neuhauser et al. 
2000; Teixeira et al. 2000; Bertout et al. 1999; van den Ancker 
et al. 1997; Kenyon & Hartmann 1995; Neuhauser et al. 1995; The 
et al. 1994) and/ or fulfilling the lithium criterion from Neuhauser 
(1997) (for individual stars see also Konig 2003; Cutispoto et al. 
2002; Neuhauser et al. 2002; Song et al. 2002; Lawson et al. 2001; 
Montes et al. 2001; Zuckerman et al. 2001; De Medeiros et al. 
2000; Gaidos et al. 2000; Strassmeier et al. 2000; Teixeira et al. 
2000; Torres et al. 2000; Zuckerman & Webb 2000; Webb et al. 
1999; Gaidos 1998; Neuhauser & Brandner 1998; Soderblom et al. 
1998; Favata et al. 1997; Micela et al. 1997; Taghaferri et al. 1997; 
Jeffries 1995; White et al. 1989). 

For the 236 pre-main sequence stars in the sample we used pre- 
main sequence evolutionary models from D'Antona & Mazzitelli 
(1994, 1997), Bernasconi (1996), Baraffe et al. (1998), Palla & 
Stabler (1999), Siess et al. (2000), Maeder & Behrend (2002) and 
Marques et al. (2008)^ to estimate their masses and ages. 
For main sequence and post-main sequence stars we used evolu- 
tionary models starting from the zero-age main sequence (ZAMS) 
from Schaller et al. (1992), Claret (2004), Pietrinferni et al. 
(2004)'*, Bertelli et al. (2008, 2009) and Marques et al. (2008). For 
stars on or above the ZAMS all models yield masses and ages based 
on luminosity and effective temperature. Stars which lie below the 
model ZAMS (mainly due to large uncertainties in the parallax, 
thus in the luminosity) are shifted towards the ZAMS, hence will 
be treated as ZAMS stars. For early-type stars this has no effect on 
the age selection (see equation 1) because they are younger than 
50 Myr even if the actual position in the HRD would be above the 
ZAMS. For most mid-F to mid-G stars age determination is dif- 
ficult and the error on the age is larger than the age value itself 
(expressing the time the star would evolve without moving in the 
HRD), thus they will not enter our list of young stars due to our 
age criterion (see equation 1). Even later late-type stars (later than 
mid-G) are already older than 50 Myr on the ZAMS. 
We used Solar metallicity for all 101628 stars in our sample with 
known magnitudes and spectral types; we obtained masses and ages 
from luminosities and temperatures. For stars with unknown lumi- 



^ Since we are working with individual stars we do not correct for statisti- 
cal biases (Smith & Eichhorn 1996). In any case, for the stars in our sample 
the Smith-Eichhorn corrected parallaxes do not differ significantly from the 
measured ones. 

■^ For some stars with unknown luminosity class we assume luminosity 
class V. For the effective temperature T^ff differences between different 



luminosity classes are modest and the error of the luminosity L is mainly 

caused by the error on the parallax. 

■^ http://www.astro.up.pt/corot/models/cesam 

■* http://albione.oa-teramo.inaf.it/ 
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nosity class, we adopt luminosity class V'. For some stars where 
also the spectral type is uncertain to ±5 sub-types, we calculate 
masses and ages for all spectral sub-types (e.g. for a G star: GO to 
G9) to derive the mean and the standard deviation of mass and age. 
We define a star to be "young" if its age is ^ 50 Myr. This limit 
is set for the following reason arising from the BSS formation sce- 
nario of runaway stars (section 1): It is desirable to identify the 
(now isolated) neutron star which was formed in the SN that re- 
leased the runaway star. This would also yield the runaway's ori- 
gin. For neutron stars the radial velocity is unknown, thus must be 
treated with a probability distribution (Tetzlaff et al. 2010). For that 
reason the error cone of the spatial motion is large and the position 
of a neutron star can be determined reliably only for a few million 
years, optimistically ~ 5 Myr. This is the maximum runaway time 
(the kinematic age) of the runaway star (as well as the neutron star) 
such that the neutron star could be identified. The latest spectral 
type on the main-sequence for stars to explode in a SN and even- 
tually become a neutron star is B3. These stars live about 35 Myr 
before they end their lives in SNe. We allow for an uncertainty of 
10 Myr and find the maximum star age (not the kinematic age) of 
the BSS runaway for which the associated neutron star should still 
be identifiable to be ~ 50 Myr. Despite that, a larger age would 
mean a longer timespan to trace back the star to identify its origin 
(if it is a DES runaway). This would cause large error bars on the 
past position of the runaway star that would make it less reliable to 
find the origin. 

However, star ages often suffer from large uncertainties due to large 
errors in distances and strong uncertainties in evolutionary models. 
That is why we choose the following criterion on the star age r* 
(the median from all evolutionary models applied) for a star to be 
young within our definition: 



both T^, + Or,* < 100 Myr and r* 



s; 50 Myr, (1) 



where cr^^* is the median deviation of t^ for different models. With 
this criterion we allow for an error of r* that is of the order of r* 
itself but also exclude stars with accurately known ages above our 
limit of 50 Myr. Unfortunately, ages of supergiants are rather un- 
certain and we miss many of them applying the criterion. For that 
reason we add all stars of luminosity classes I and II as well as stars 
of luminosity class III earlier than AO because they could certainly 
be younger than 50 Myr. Moreover, the classical definition of run- 
away stars by Blaauw (1961) includes stars up to B5 of luminos- 
ity classes IV and V. These are added as well. Finally, our sample 
contains 7663 young stars. Their Hipparcos numbers and common 
names as well as ages, masses and spectral types are listed in Ta- 
ble 1 (the full table will be available via the VizieR database soon 
after publication). Since 2250 stars in the table enter our list only 
due to their spectral type and luminosity class, we list only their 
spectral type and luminosity class as given in the literature. 



Table 1. Ages r* (in Myr), masses Af* (in solar masses) and spectral types 
for 7663 young stars (sorted by their HIP number), r* and Af* are medi- 
ans obtained from different evolutionary models (see section 2.1). For 2250 
stars only the spectral type is given since models infer a larger age; however 
they are probably also young (as inferred from spectral type and luminosity 
class). Here the first five entries are shown, the full table will be available 
via VizieR soon after publication. 



HIP 


other name 


r* [Myr] 


Af* [Mq] 


SpT 


32 


HD 224756 


10.0 ± 7.6 


2.8± 0.1 


B8 


85 


CD-25 16747 


39.8 ± 28.8 


1.8± 0.1 


A2 


106 


HD 224870 


- 


- 


G7II-III 


124 


HD 224893 


44.7 ± 5.4 


7.6 ± 0.4 


FOIII 


135 


HD 224908 


45.7± 22.9 


l.Oit 0.1 


G5 



the kinematic centre of the stars in our sample as follows. We cal- 
culate spatial velocity components of the 4195 stars with complete 
kinematic data (in a right-handed coordinate system with the x axis 
pointing towards the Galactic centre and y is positive in the direc- 
tion of galactic rotation) corrected for Galactic differential rotation 
using Keplerian orbits: 



U= U^Urot, 

V= V -Vrot + VQ^rot, 

w = w. 



(2) 



where U , V and W are the heliocentric velocity components (in 
the X, y and z direction, respectively) and Urot and Vrot the com- 
ponents of the rotational velocity of the star moving around the 
Galactic centre. VQ^rot = 225kni/s is the rotational velocity of 
the Sun around the Galactic centre. To avoid significant contamina- 
tion of high velocity stars, we exclude stars with v > 50 km/s (that 
is approximately two times the median of {) = v f/^ + 1/^ + W'^). 
Fitting a Gaussian to each velocity component, we find 

(Uq,V0,Wq) = (10.4 ±0.4, 11.6 ±0.2, 6.1 ±0.2) km/s. 

(3) 
The agreement of our LSR with the classical value of Delhaye 
(1965) [(9,11,6) km/s] is remarkable. In comparison with the most 
recent value published by Aumer & Binney (2009) [(9.96 ± 
0.33,5.25 ± 0.54,7.07 ± 0.34) km/s], the V component differs sig- 
nificantly; however, Aumer & Binney (2009) obtained their results 
by examining the correlation between LSR and colour B — V, i.e. 
star age, and extrapolating the curves to zero velocity dispersion 
and also ignoring stars with B — V^O (see also Dehnen & Binney 
1998) because they are probably not yet mixed. Such stars are over- 
abundant in our sample of young stars. From Figure 3 in Dehnen & 
Binney (1998) one may easily see that our result agrees well with 
their findings. 

In the following section we will use equation 3 to correct the veloc- 
ities for Solar motion. The peculiar velocity is the velocity of a star 
corrected for Solar motion and Galactic rotation. 



2.2 Kinematics 

The solar motion with respect to a specific star sample, i.e. the Lo- 
cal Standard of Rest (LSR), depends upon the age of the stars in the 
sample (e.g. Mihalas & Binney 1981). For that reason, we derive 



^ For the effective temperature T^ff differences between different lumi- 
nosity classes are modest and the error of the luminosity L is mainly caused 
by the error on the parallax. 



3 YOUNG RUNAWAY STARS 

Classically, a star is assigned to be a runaway star if its peculiar 
space velocity Vp^c exceeds 40 krn/s (Blaauw 1961). Later on, Vit- 
richenko et al. (1965) and Cruz-Gonzalez et al. (1974) used ra- 
dial velocities Vr alone to investigate the population of runaway 
stars defining a more general definition of |wr,pec| > Vcrit,iD with 
VcritAD ~ 3(T, whcrc a is the mean velocity dispersion in one di- 
mension of the low velocity stars. Based on this definition Moffat 
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et al. (1998) investigated tangential velocities vt to identify run- 
away stars. We now want to combine and extend all these methods. 

3.1 Runaway stars identified from their peculiar space 
velocity 

Runaway stars were firstly described as the stars that are responsi- 
ble for the longer tail in the velocity distribution such that it is not 
sufficiently describable with one Maxwellian distribution (Blaauw 
1961). Stone (1979) generalised the definition such that runaway 
stars are the members of the so-called high velocity group. These 
are stars with large peculiar velocities that can be represented by an 
additional Maxwellian distribution. The other Maxwellian distribu- 
tion incorporates stars with lower velocities, thus the low velocity 
group (normal Population I stars). He pointed out that by applying 
a velocity cutoff to identify the runaway stars, a certain fraction of 
them would be missed. 

However, this issue can only be handled for the determination of 
space frequencies (see Stone 1991) and a velocity cutoff is still in- 
evitable for the identification of runaway star candidates. To obtain 
a reasonable cutoff, we fit the distribution of the peculiar space ve- 
locities Vpec of the sample stars (4180 with full kinematics) with 
two Maxwellians (Fig. 1). We evaluate the velocity errors utilising 
a Monte-Carlo simulation varying tt, fi^, fJ-s and Vr within their 
confidence intervals.'' We find that 

(TL =9.2±0.2km/s, 

an =24.4±1.2km/s, (4) 

fn =27.7 ±1.9%, 

where a^ and an are the average velocity dispersions of the low 
and high velocity groups, respectively, and fn is the relative fre- 
quency of the high velocity group. The derived dispersions are 
in agreement with those found by Stone (1979) (with an being 
slightly smaller) whereas fn is smaller; however published values 
for /h vary from 34 ± 14 % (Blaauw 1961, corrected - see Stone 
1991) to 55 ± 12 % (Stone 1979). Furthermore, Stone's star sam- 
ple contains a much smaller number of stars than ours. In addition, 
to make sure that the low mass stars in our sample do not distort 
the results, we check whether the outcome differs from a subsam- 
ple comprising only O and B type stars as well as Wolf-Rayet stars 
(2368 stars with full kinematics). Since we do not find significant 
differences we conclude that young stars, no matter if of low or high 
mass, share the same kinematic properties. The Maxwellian func- 
tions of both groups intersect at 28kin/s. Following Stone (1979), 
a star is a probable member of the high velocity group if 



> 28 km/s. 



(5) 



For that reason this will be our velocity cutoff defining runaway 
stars. In theory, with this definition we are able to identify 73 per 
cent of the high velocity group members while the contamination 
of low velocity stars is nine per cent. 

We perform a Monte-Carlo simulation varying the observables vr, 
Ha, Mi5 3nd Vr within their uncertainty intervals and evaluate the 
probability of a star being a runaway star (the probabilities are 
given in Table 4). For 972 stars, the probability is higher than 50 
per cent. Allowing for a nine per cent contamination of low veloc- 
ity stars, this means a probability criterion of 50 per cent allows us 
to identify 78 per cent of the high velocity members. 

^ Note that the velocity errors are not symmetric due to the inverse propor- 
tionahty regarding tt. 




40 60 

Vpec [km/s] 



100 



Figure 1. Distribution of the peculiar 3D space velocity Vp^c (shaded his- 
togram). The dashed curve shows the distribution for the low velocity group 
whereas the dashed-dotted curve is for the high velocity group. The two 
curves intersect at Vpec = 28 kin/s. The total distribution as the sum of the 
two is represented by the full Hne. Considering that different moving groups 
and associations of stars are present, thus both, the low and high velocity 
groups of stars, might be actually better represented by a superposition of 
many Maxwellians with slightly different parameters, the fit iXred ~ 1-57) 
is satisfying. 



3.2 Runaway stars identified from U, V, W, tlieir radial and 
tangential velocities or proper motions 

In addition to the peculiar 3D space velocities Vp^c we investigate 
their ID components U, V and W separately to identify poten- 
tially slower high velocity group members which may show an 
exceptionally high velocity in only one direction. For the same 
reason we also investigate the peculiar radial velocities Vr,pec- 
45 per cent of the stars in our sample do not have radial velocity 
measurements available. Among these cases, the only way to 
identify runaway candidates is to use their peculiar 2D tangential 
Galactic velocities vt,pec or its ID components which are the 
peculiar proper motion in Galactic longitude ^i^pec and Galactic 
latitude iib,pec- To make the velocities comparable, we transfer the 
proper motions into ID velocities «i,pec = 4.74 • pii,pec/7r and 

Vb,pec = 4.74 ■ fJ,b,pec/T'- 

All velocity distributions contain the two velocity groups 
of stars (section 3.1) and can be fitted with bimodal functions 
(Gaussians for the ID cases U, V, W, Vr,pec, vi.pec and Vb.pec, 2D 
Maxwellians for the 2D case vt,pec)- Table 2 lists the fitting results 
adopting fn ~ 27.7 ± 1.9 % (see section 3.1). 
The velocity dispersions of the high velocity group are consistent 
with an isotropic velocity distribution arising from the runaway 
producing mechanisms (see appendix A). Moreover, the low 
velocity group dispersions are in good agreement with those of 
young disk stars (e.g. Delhaye 1965; Mihalas & Binney 1981). 
Since the velocity distribution of the low velocity group is not 
isotropic the one of the high velocity group cannot be isotropic 
either (see appendix A). Thus, we cannot simply translate the cri- 
terion given by equation 5 into the ID case (\X\ > 28/\/3km/s, 



where X = U,V, W, Vr 



,VU; 



or «(,,pec). In addition, such a 



ID criterion would lead to a contamination of low velocity group 
stars among the identifications of approximately 50 per cent e.g. in 
the U and V components. 
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Table 2. Fit results and curve intersection points for the different velocity components. The first Column gives the velocity component examined. Columns 2 
and 4 give the centre velocities Dd, and v^-h (for the ID cases, i.e. Gaussian fit) and Columns 3 and 5 give the velocity dispersions crj^ and ctjj of the low 
and high velocity groups, respectively. All errors are formal la errors. In Column 6 the intersection points of the curves representing the low and high velocity 
groups are given (for the ID cases two intersection points exist (negative and positive sides of the distribution), they are approximate since the distribution is 
not exactly symmetric). The last Column names the Figure in which the particular distribution is shown. 





VcL [km/s] 


CTi [km/s] 


VcH [km/s] 


dH [km/s] 


intersection [km/s] 


Fig. 


u 


-0.3± 0.3 


10.7± 0.3 


4.1 ± 1.5 


26.3 ± 1.0 


±23 


2(a) 


V 


0.0 ± 0.2 


10.7± 0.1 


-4.7± 1.0 


24.2 ± 1.0 


±23 


2(b) 


w 


O.Oit 0.1 


5.3± 0.1 


-3.8± 0.5 


17.3 ± 0.6 


±12 


2(c) 


'^r,pec 


0.2± 0.3 


11.9± 0.2 


-4.9± 1.1 


28.5 ± 0.9 


±25 


3 


Vt,pec 


- 


7.4± 0.1 


- 


21. 7± 0.6 


20 


4 


'"l,pec 


-2.7± 0.1 


8.9± 0.1 


-2.6± 0.7 


27.3 ± 0.8 


±19 


5(a) 


'"b,pec 


0.6± 0.1 


5.3± 0.1 


-3.3 ± 0.4 


18. 3± 0.7 


±11 


5(b) 





100 -100 



[' [km/s 

(a) U, xl^a = 0-73 



V [km/s 



W [km/s] 

W W, x?ed = 0-53 



Figure 2. Distribution of U, V and W (shaded histograms) - the ID components of Vpsc- The dashed curves represent the low velocity group while the 
dashed-dotted curves display the high velocity group. The two curves intersect at (7 ft: ±23 km/s, V ~ ±23 km/s and VF ~ ±12 km/s, respectively. The 
total distribution as the sum of the two is represented by the full line. 



140 




-50 

Vr,pec [km/s] 



Figure 3. Distribution of the peculiar ID radial velocity Vr 



100 



(shaded 



histogram). The dashed curve represents the low velocity group while the 
dashed-dotted curve displays the high velocity group. The two curves inter- 
sect at Vr^pec ~ ±25 km/s. The total distribution as the sum of the two is 
represented by the full line (Xred ~ 0-33). 



For those reasons, as for the peculiar spatial velocity Vpec, we 
choose the intersection points of the curves also for the ID cases 
as well as Vt,pec to define the runaway criteria. In Table 3, we 



specify the selection criteria, theoretical expectations concerning 
the possible identifications as well as the number of runaway 
star candidates identified (stars with a runaway probability 
higher than 50 per cent). Allowing for the specific contamination 
of low velocity group members, the number of identifications 
is generally in agreement with our theoretical predictions (Table 3). 

With 972 runaway star candidates found in section 3.1 and 
1381 identifications in Table 3 (Nnew), we find a total of 2353 run- 
away star candidates with a runaway probability higher than 50 per 
cent regarding at least one velocity investigated (Table 4, the full 
table will be available via VizieR soon after publication). 



3.3 Stars with higher pecuUar velocities compared to their 
neighbourhood or surrounding OB association/ cluster 

Still, some high velocity group stars may not yet have been 
identified. For that reason we look for additional stars which show 
a different motion compared to their neighbouring stars. Since stars 
in clusters and associations share a common motion, runaway stars, 
i.e. stars that experienced some interaction (see section 1), can be 
identified through deviations from the common motion, especially 
if the velocity vector points towards a different direction than the 
cluster mean motion. From the previous investigations we can be 
sure that we identified all the runaway stars with high peculiar 
velocities. The most important criterion now is the direction of a 
star's velocity vector compared to its neighbouring stars. 
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Table 4. Runaway probabilities for 2353 runaway star candidates as found in the previous section. Columns 3 to 10 list the individual probabilities P for each 
velocity component. We regard stars with P ^ 0.50 in at least one velocity as runaways. The peculiar space velocities ?^pec and peculiar tangential velocities 
iJt.pec are given in the last two Columns. Here the first five entries are shown, the full table will be available via VizieR soon after publication. 



HIP 



other name 



Pv„ 



Pv, 



"t.i 



Pn^; 



Pv. 






[km/s] vt,pec [km/s] 



85 
135 
145 
174 
278 



CD-25 16747 
HD 224908 

29Psc 
HD 240475 
HD 225095 



0.89 

0.13 
1.00 
0.45 



0.93 

0.00 
0.52 
0.13 



0.47 

0.00 
1.00 
0.23 



0.22 

0.95 
0.00 
0.04 



0.01 

0.31 
1.00 
0.51 



0.96 
0.99 
0.00 
0.04 
0.03 



0.52 
0.00 
0.00 
0.04 
0.02 



0.99 
1.00 
0.00 
0.00 
0.03 



52.91 

22.4 
59.6 
27.3 



h22.8 
-35.2 


c:9 7+24.8 




22.1+1-? 


+3.9 
-4.1 


1 9+^° 


+4.8 
-5.2 


6.9t?:°o 


+8.9 
-7.1 


9.5t?:« 




40 60 

Vt,pec. [km/s] 



100 



Figure 4. Distribution of the peculiar 2D tangential Galactic velocity vt,pec 
(shaded histogram). The dashed curve represents the low velocity group 
while the dashed-dotted curve displays the high velocity group. The two 
curves intersect at Vr.pec = 20 km/s. The total distribution as the sum 
of the two is represented by the full line (x^ed ~ 2.58, see comment to 
Fig. 1). 



Table 3. Runaway selection criteria |i;| > Vcrit for the different velocity 
components. The limits coiTespond to the intersection points of the curves 
representing the low and high velocity groups. Columns 3 and 4 give the 
theoretical expectations concerning possible identifications of high velocity 
group members (the fraction of high velocity group stars /^^ tj, satisfying 
the selection criteria) and the contamination of low velocity group mem- 
bers (the fraction of low velocity group stars fc^th that also satisfy the se- 
lection criteria) while in the last two Columns the number of runaway star 
candidates A'^ and the number of new identifications Nnew (compared to 
section 3.1 and previous lines in the table) is quoted. Note that for U, V, 
W and Vr,pec only 4180 stars could be analysed whereas for vt^pec, iJi^pec 
and vi,p^i. the whole sample of 7663 stars was used. 



Vcrit [km/s] fia,th [%] fc,th [%] 



N 



Nn 



\u\ 
\v\ 

\W\ 

|'^r,pec| 

\'^l,pe.c 
\'^b,pec 



23 
23 
12 
25 
20 
19 
11 



39 
35 
50 
39 
66 
49 
56 



17 


776 


70 


19 


452 


43 


12 


609 


190 


20 


588 


12 


9 


1513 


768=' 


18 


1162 


33'' 


15 


1266 


265<= 



^ 643 of them without Vr measurements 
'' 33 of them without v,- measurements 
'^ 265 of them without Vr measurements 




(a) fi.pec, Xred = 0-47 



(b) Vt^pec, X 



2 
red 



1.06 



Figure 5. Distribution of the ID components of the peculiar 2D tangen- 
tial Galactic velocity Vt^pec- *'!,pec and f^^pec (shaded histograms). The 
dashed curves represent the low velocity group while the dashed-dotted 
curves display the high velocity group. The two curves intersect at vip^c ~ 
±19 km/s and vij^^c ~ ±11 km/s, respectively. The total distribution as 
the sum of the two is represented by the full line. 



We select a sphere with a diameter of 24 pc^ around each individ- 
ual star to define its neighbourhood. All sample stars within this 



sphere are chosen as comparison stars. We calculate the vectors 
of the peculiar velocities Vp^c = {U, V, W) and the peculiar tan- 
gential velocities Vt,pec = (^;,pec, ^6, pec) varying the observables 
within their confidence intervals. The neighbourhood velocity 
Vneigh IS defined as the median velocity of the comparison stars. 
We define the runaway criterion such that w* must not lie within 
the 3(7 error cone of v„eigh (Fig. 6). Varying the observables 
within their confidence intervals, we obtain runaway probabiUties 
for each star (10000 Monte-Carlo runs). 

Examining Vp^c and vt,pec, we find no stars with a probability 
of at least 50 per cent for being a runaway star under the above 
definition. 

In addition, we compare the velocity vectors of each OB as- 
sociation and cluster listed in Tetzlaff et al. (2010)*^ with the ve- 
locity vectors of each individual star situated within the association 
boundaries as listed in Tetzlaff et al. (2010) (assuming the associa- 
tions/ clusters are spherical, see Fig. 7 for the Orion OB 1 associa- 
tion as an example). From Vp^c, we identify 126 additional runaway 
star candidates. Another 58 runaway star candidates are identified 
from their Vt,pec- All runaway star candidates (192 for Vp^c, 221 



^ This is twice the median extension of all associations listed by Tetzlaff 
et al. (2010). 



** We exclude Cas-Tau from the analysis owing to its large size (cf. de 
Zeeuw et al. 1999, their Table 2). 
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Table 5. Runaway probabilities for runaway star candidates found by comparison with OB associations and clusters (as listed in Tetzlaff et al. 2010). 192 
runaway stars are identified from Vp^c and 221 from vt.pec- 124 stars are included in both lists. Columns 3 and 4 list the runaway probability P as well as 
the association/ cluster. The absolute velocity values are given in columns 5 and 6. Errors correspond to 68 per cent confidence. The last Column indicates 
whether the star was already identified as a runaway star in sections 3.1 and 3.2 ("prev" for previous identification; "new" for new identification). Here the first 



entries for Vp^c 


and Vt^pec, respectively, are shown, the full table will be available 


via VizieR soon 


after publication. 






HIP 


other name 


P 


Assoc./ cluster 


Vpec [km/s] 


vt,pec [km/s] 


new ident. 










from Vpec 










490 


HD105 


1.00 
0.99 


13 Pic-Cap 
ABDor 


ll.oi?;? 


10.6±?;? 


new 




1481 


HD 1466 


1.00 
1.00 
1.00 


Tuc-Hor 

13 Pic-Cap 

ABDor 


10.81?;? 


10.71?;? 


new 




1623 


HD 1686 


1.00 
1.00 


Tuc-Hor 
P Pic-Cap 


27.7t\l 


27.5t?;? 


prev 




1803 


BECet 


1.00 
1.00 
1.00 


Tuc-Hor 

13 Pic-Cap 

ABDor 


27.41?;? 


27.1±?;? 


prev 




1910 


GSC 08841-00065 


0.65 


13 Pic-Cap 


13.1±^;? 


IS.Ol^;? 


new 










from vt,pec 










135 


HD 224908 


0.96 


ABDor 


- 


22.ltl-? 


prev 




439 


HD 225213 


1.00 
1.00 


Tuc-Hor 
ABDor 


iii.3±?;^5 


109.51?;^ 


prev 




490 


HD 105 


1.00 
1.00 


ABDor 
13 Pic-Cap 


ll.ol?;? 


I0.6t?;? 


new 




544 


V 439 And 


1.00 
1.00 


Tuc-Hor 
P Pic-Cap 


ll.9ii;? 


ll.ll?;? 


new 




1128 


HD967 


1.00 
1.00 


Tuc-Hor 
13 Pic-Cap 


9.2.2tll 


77.s+J;l 


prev 




Figure 6. Definition for identifying runaway stars by comparing with neigh- 
bouring stars. The neighbourhood is defined as a sphere around an individ- 
ual star with a diameter of 24 pc (footnote 7) or as the OB association/ 
cluster (assumed to be spherical) inside which the star currently lies. The 
neighbourhood velocity v„eigh is given by the median velocity of the stars 
within the sphere or the mean motion of the association/ cluster member 
stars as listed in Tetzlaff et al. (2010), respectively. If the velocity vector 
Vi, points into the grey shaded region, the star clearly moves away from its 
neighbouring stars and is hence a runaway star. The grey shaded area lies 
outside the 3cr error cone of "CVieigh (dotted lines mark la). Note that the 
length of Vir and 'O'neigh (Vpec or vt,pec) is not important here. 



for Vt,pec\ 124 included in both) are listed in Table 5 (the full table 
will be available via VizieR soon after publication). 



0) -15- 




-25 



205 210 

Galactic longitude I [°; 



Figure 7. Motion of Ori OBI member stars. The red cross marks the cen- 
tre of the association, the dotted ellipse its boundaries (assuming spheri- 
cal shape). The red thick arrow shows the mean peculiar tangential motion 
of the whole association (Blaha & Humphreys 1989; Brown et al. 1999; 
Dambis et al. 2001; see also Tetzlaff et al. 2010), green squares are run- 
away star candidates satisfying the criterion defined by Fig. 6 whereas blue 
stars mark runaway star candidates already defined by their large peculiar 
velocity (section 3). The length of the aiTows is scaled with distance to in- 
dicate tangential velocities. 
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Table 6. Additional young stars situated well outside any OB association/ cluster and the Galactic plane (see section 3.4 for criterion), i.e. runaway star 
candidates. Columns 3 to 4 give the distance z to the Galactic plane as well as the star's velocity Wpec in this direction. The absolute velocity values are given 
in columns 5 and 6. Errors correspond to 68 per cent confidence. Although some stars would be consistent with z = within 2cr and might be outliers we 
include them as potential runaway stars to not miss them only due to their large error bars. 



HIP 


other name 


z[pc] 


5805 


HD 7598 


.77+280 
*' '-350 


11242 


HD 14920 


-6091^10 


50684 


RSSex 


t;c:c: + 135 
000-245 


54769 


HD 97443 


5^7+111 


56473 


90 Leo 


/IOO+320 
*°°-280 


70000 


HD 125504 


t-44+300 



Wpec [km/s] 



[km/s] vt,pec [km/s] 



-16.3 



-3.2"' 



-10.0 
-8.0 



12.0 



-1-5.9 
10.1 



-10.2 



10.0 

fS.O 



14.3 
25.4 



+ 7.9 
10.1 
-1-6.9 
9.1 



13.1 
15.9 

6.6 

5.7 

8.0 

20.0 



-1-12.9 

13.1 

-F9.9 

16.1 

-1-5.0 

7.0 

-1-4.0 

6.0 

-f6.0 

8.0 

-1-3.9 

1 



3.4 Stars outside OB associations/ clusters and tlie Galactic 
plane 

As runaway stars were ejected from their birth site, i.e. its host 
OB association/ cluster or the Galactic plane, they are supposed to 
be isolated (outside any OB association/ cluster and the Galactic 
plane). Thus, we look for young stars that are clearly outside any 
OB association/ cluster listed in Tetzlaff et al. (2010) (outside three 
times the association radius which corresponds to approximately 
3a) and probably outside the Galactic plane (z > 500 pc®). 
72 stars are situated well outside any OB association/ cluster and 
the Galactic plane; six of them were not identified as runaway star 
candidates in the previous sections. Those six are listed in Table 6. 
If they did not form in isolation they are runaway stars. 



3.5 Comparison with other authors 

We compare our sample of runaway star candidates with lists from 
other authors (most important sources: Blaauw 1961; Bekenstein 
& Bowers 1974; Cruz-Gonzalez et al. 1974; Stone 1979; Gies & 
Bolton 1986; Leonard & Duncan 1990; Conlon et al. 1990; Philp 
et al. 1996; Moffat et al. 1999; Hoogerwerf et al. 2001; Mdzinar- 
ishvili & Chargeishvili 2005; de Wit et al. 2005; Martin 2006). 24 
proposed runaway candidates which satisfy our initial sample cri- 
teria (HIP star, TT — CT^ ^ 1/3 mas, equation 1), i.e. are contained 
in our young star sample, were not identified as runaway stars by 
us. They are listed in Table 7 along with the respective publication 
source. 

We add three of the missing stars, namely HIP 3881 (=35 And), 
HIP 48943 (=0Y Hya) and HIP 102195 (=V4568 Cyg) as well as 
HIP 26241 (=(, Ori) due to its DES origin (see individual discussion 
in appendix B). 

As indicated by the discussion in appendix B, several problems may 
lead to mis- or non-identification of runaway stars. The major issue 
here is certainly the distance of a star that highly affects the calcu- 
lated velocities. Since we used precise Hipparcos parallaxes (while 
previous studies often used ground-based distances) our results are 
not significantly influenced by that. Moreover, we derive runaway 
probabilities accounting for the errors on all observables instead of 
evaluating only a single velocity. 

Another problem arises from the multiplicity of stars, e.g. 390 stars 
in our young star sample are spectroscopic binaries (Pourbaix et al. 



^ This number is derived from twice the low velocity dispersion in the z di- 
rection (ai 10 km/s) and an age of 50 Myr (our age limit, see section 2.1). 
We do not use the individual ages of the stars owing to the large uncertainty 
inferred from different evolutionary models. 



Table 7. Runaway star candidates from Hipparcos proposed in the literature 
that are absent from our sample of 2543 runaway stars. 



HIP 



Ref. 



102195 Blaauw (1961); Bekenstein & Bowers (1974) 
35412 Cruz-Gonzalez et al. (1974) 

40328 Cruz-Gonzalez et al. (1974) 

81736 Cruz-Gonzalez et al. (1974) 

18246 Stone (1979) 

100214 Stone (1979) 

114104 Stone (1979) 

117221 Stone (1979) 

67279 Leonard & Duncan ( 1 990) 

3881 Hoogerwerf et al. (2001) 

20330 Hoogerwerf et al. (2001) 

38455 Hoogerwerf et al. (2001) 

48943 Hoogerwerf et al. (200 1 ) 

86768 Hoogerwerf et al. (200 1 ) 

92609 Hoogerwerf et al. (200 1 ) 

103206 Hoogerwerf et al. (2001) 

69640 Mdzinarishvili & Chargeishvili (2005) 

73720 Mdzinarishvili & Chargeishvili (2005) 

80338 Mdzinarishvili & Chargeishvili (2005) 

90074 Mdzinarishvili & Chargeishvili (2005) 

104361 Mdzinarishvili & Chargeishvili (2005) 

1415 de Wit et al (2005) 

37074 de Wit et al. (2005) 

97757 de Wit et al. (2005) 

2009; Pedoussaut et al. 1988). We obtained radial velocities from 
catalogues which typically list the systemic radial velocity. 



4 SUMMARY AND CONCLUSIONS 

We analysed the distributions of the peculiar velocities of 7663 
young Hipparcos stars (4180 with full kinematics) in three, two 
as well as one dimension to identify members of the high ve- 
locity group, i.e. stars which show different kinematics than nor- 
mal Population I objects and hence experienced some interaction 
(BSS or DES) which provides them with an additional velocity. 
As expected, the velocity component due to runaway formation is 
isotropic. 

Performing Monte-Carlo simulations varying the observables n, 
Ha > f^s and Vr within their uncertainty intervals we assigned each 
star a probability for being a member of the high velocity group, 
i.e. a runaway star. We did that for the 3D velocity Vp^c, its ID 
components U, V and W, the ID radial velocity Vr,pec as well as 
the 2D tangential velocity vt,pec and its ID components vi^p^c and 
Vb,pec to identify as many members of the high velocity group as 



A catalogue of young runaway stars 9 




Figure 8. Distributions ofvt^pec and Vpec- Light grey liistograms represent 
tlie wliole star sample whereas dark grey histograms show the distributions 
of runaway star candidates. The velocity distributions of the high velocity 
group is drawn with dashed-dotted lines (cf. Figs. 4 and 1). 



possible, also those with a relatively low velocity. We found 2353 
runaway star candidates, i.e. stars for which the runaway probabil- 
ity was higher than 50 per cent in at least one velocity component 
examined. The contamination of low velocity members is about 20 
per cent at most. 

However, those high velocity group members with very small ve- 
locities (as inferred for a Maxwellian distribution they exist) could 
still not be identified due to the velocity thresholds set. Therefore, 
we compared the velocity vector (in 3D, Vpec, and 2D, vt.pec) of 
each individual star with the one defined by the stars in its neigh- 
bourhood or its surrounding OB association/ cluster (Vneigh)- If the 
velocity of the star clearly pointed away from Vnetgh (see Fig. 6), 
the star was assigned being a runaway star. 1 84 new identifications 
were made. 

Additionally, we find six additional young stars which are situated 
well outside any OB association/ cluster and the Galactic plane. 
Finally, we compared our list of runaway star candidates with pre- 
viously published lists and added four stars (see discussion in ap- 
pendix B). 

hat gives us a total of 2547 runaway star candidates (with a con- 
tamination of low velocity group members, i.e. normal Population 
I stars, of 20 per cent at most). Thus the runaway frequency among 
young stars is approximately 27 per cent, in agreement with our 
theoretical expectations. 

Fig. 8 shows the distribution of the peculiar space velocity Vpec 
and the peculiar tangential velocity vt^pec (as Figs. 1 and 4) with 
the subsample of runaway candidates in dark grey. The number of 
runaway star candidates is somewhat higher than the number of 
stars belonging to the high velocity group, especially in the range 
of medium velocities where the two curves intersect, as to be ex- 
pected. However, with this somehow conservative selection, we can 
be sure that we did not miss actual runaway stars. Using our com- 
bined selection we also identified runaway star candidates with rel- 
atively low velocities which certainly exist but would have not been 
identified by investigating only one velocity component or the ab- 
solute velocity. 

We will provide Tables 1, 4, 5 and 6 within an electronic cat- 
alogue available via VizieR. 
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APPENDIX A: THE VELOCITY DISPERSION OF THE 
HIGH VELOCITY GROUP OF STARS 

For normal Population I stars, i.e. members of the low velocity 
group, the velocity dispersion in the z direction, aw, is smaller 
than the velocity dispersion in the x and y directions (in the Galac- 
tic plane) due to the Galactic potential that attracts the stars onto 
the Galactic plane. Since all stars initially belong to the low veloc- 
ity group there must be a difference between au/v and aw also 
for high velocity group members. Since the velocity distribution in 
each direction is Gaussian, 



2 


2 , 2 
= '^L,U + <^x, 


^H,V = 


= (^L,V+(^x, 


'^H,W = 


= (^L,W + (^x 



(Al) 



where ax is the velocity dispersion due to runaway formation. Ini- 
tially, we assume the additional velocity to be isotropic. For the low 
velocity group we find that 



aL,w 
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Thus, 
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<^H,w = <^H,u/v ~ 10 km/s • (tl, IV — 25 km /s . 
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With au^u/Y ~ 25 km/s and aL,w ~ 5 km/s, it follows that 
aH,w ~ 23 km/s theoretically. However, since the runaway for- 
mation occurred some time in the past (for BSS runaways in the 
sample this timespan might be comparable with the age of the star 
before the supernova), the Galactic potential makes an impact on 
the higher velocities. For that reason, we expect to measure a some- 
what lower value aw,H for the high velocity group than the pre- 
dicted value of 23 km/s, thus the found value of aw,H ~ 17 km/s 
is in good agreement with our predictions and we conclude that 
runaway formation leads to an additional velocity that is isotropic. 



APPENDIX B: AN INDIVIDUAL DISCUSSION ON 
RUNAWAY STARS FOUND IN THE LITERATURE 

The classical runaway HIP 102195 was not identified by us since 
its peculiar spatial velocity is rather small {Vpec = 16. GI^*^ km/s). 
Note that Blaauw (1961) also quote a small velocity (~ 23 km/s). 
As proposed by the authors, HIP 102195 apparently originated 
from the Lacerta OBI association. However, using 3D data, we 
cannot confirm this origin but instead find that in 10.0 % of 10000 



Monte-Carlo runs the star's position is located within the bound- 
aries of Cygnus OB7 about ll.llgf Myr in the past which is in 
excellent agreement with the association age of 13 Myr (Uyaniker 
et al. 2001). Due to the large number of parameters involved (po- 
sition and velocity of the star and association) the fraction of suc- 
cessful runs is expected to be small (cf. Hoogerwerf et al. 2001; 
Tetzlaff et al. 2010). Hence, we include this star into our runaway 
star sample. 

The twelve stars identified by Cruz-Gonzalez et al. (1974), Stone 
(1979) and Mdzinarishvili & Chargeishvili (2005) are not re- 
identified by us simply because the authors used photometric 
distances which are systematically too large, thus generating 
large peculiar velocities (this can be directly seen from compar- 
ison between Columns 5 and 6 of Table 1 in Mdzinarishvili & 
Chargeishvili 2005), whereas we used parallactic distances to de- 
termine peculiar velocities. 

HIP 67279 was included by Leonard &. Duncan (1990) owing to 
its large distance from the Galactic plane of 2 = 1 kpc (accord- 
ing to the author's definition of a runaway star z must be larger 
than 20 km/s times the main-sequence lifetime of the star, cf. foot- 
note 9). The photometric distance of 1.1 7 kpc that the authors 
adopted from Kilkenny et al. (1975) is however once more too large 
and the actual distance from the Galactic plane derived from the 
parallax (parallactic distance is 472lijQ3 pc) is z = 3971j3q pc. 
With this z and an age r* = 0.3 ± 2.1 Myr as inferred from evo- 
lutionary models (see section 2), HIP 67279 would need a vertical 
velocity component as large as, W = 1300 km/s to have origi- 
nated from the Galactic plane. Similarly, de Wit et al. (2005) iden- 
tified the three stars listed in Table 7 only from their distance to the 
Galactic plane again using photometric distances. With the better 
parallactic distances they do not satisfy the criterion applied by the 
authors (z > 250 pc). 

Four of the seven runaway candidates listed by Hoogerwerf et al. 
(2001) (HIP 20330, 86768, 92609, 103206) are not recognised as 
runaway stars by us due to different input data (especially tt) as 
Hoogerwerf et al. (2001) used the old Hipparcos reduction (Perry- 
man et al. 1997) (smaller tt in all four cases) whereas we used the 
latest published data by van Leeuwen (2007). In the case of HIP 
48943, the radial velocity adopted by Hoogerwerf et al. (2001) of 
Vr — 39.0 ± 5.0 km/s differs from our value of Vr = 29.6 ± 
3.6 km/s (Kharchenko et al. 2007) resulting in different peculiar 
space velocities (30.7^5'^ km/s and 22.6t.^'l km/s, respectively). 
Since HIP 48943 is an astrometric binary (Makarov 2007), we ac- 
count for uncertainties in the radial velocity and include it into our 
sample of runaway star candidates. For another one, HIP 38455, 
Hoogerwerf et al. (2001) adopted Vr = —31.0 ± 5.0 km/s. Ac- 
cording to Haefner & Drechsel (1986) HIP 38455 is a spectroscopic 
binary with a systemic radial velocity of 29.5 km/s. This signifi- 
cantly different radial velocity changes the peculiar spatial velocity 
dramatically (vr — —31 km/s: Vpec = 47 km/s, Vr = 29.5 km/s: 
Vpec = 14 km/s), hence the star is not a runaway. Moreover, 
Hoogerwerf et al. (2001) corrected the velocities for Solar motion 
using the LSR published by Dehnen & Binney (1998) which results 
in somewhat larger space velocities than ours (which is why we did 
not find some of their runaway stars) but does not accurately reflect 
the motion of young stars relative to the Sun (see section 2). For 
HIP 3881 Hoogerwerf et al. (2001) suggested a birth association 
(Lacerta OBI). We find that in 1.8 per cent of 10000 Monte-Carlo 
runs the star's position coincided with the boundaries of Lacerta 
OBI (6.3t'i:2 Myr in the past). Like HIP 102195 (see above), we 
include this star into our runaway star sample. In addition, we in- 
clude HIP 26241 (= t Ori, highly eccentric spectroscopic binary) 
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since it was very probably part of a former triple system (together 
with the classical runaways HIP 24575 = AE Aur and HIP 27204 
= ji Col), thus ejected via DSS from the Trapezium cluster (e.g. 
Hoogerwerf et al. 2001) and a member of the high velocity group. 



